Inelastic X-ray scattering in correlated (Mott) insulators. 
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We calculate the inelastic light scattering from X-rays, which allows the photon to transfer both 
energy and momentum to the strongly correlated charge excitations. We find that the charge transfer 
peak and the low energy peak both broaden and disperse through the Brillouin zone similar to what 
is seen in experiments in materials like Ca2Cu02Cl2. 

PACS numbers: 78.70.Ck, 72.80.Sk, 78.66.Nk, 71.30.+h, 71.27.+a 



The dynamics of electrons in strongly correlated sys- 
tems is far from well-understood. In a Mott insulator, 
correlations split a single band into a lower and an upper 
Hubbard band separated by a Mott gap. Many exper- 
imental probes have focused attention on exploring the 
detailed nature of the lower Hubbard band from which 
electrons may be excited using angle-resolved photoemis- 
sion (ARPES) for example, but the structure and sym- 
metry of the upper Hubbard band and the relaxational 
dynamics of electrons populated into it remains largely 
unexplored. 

Inelastic X-ray scattering Jl) has attempted to ad- 
dress this issue on a number of correlated insulators 
such as La 2 Cu0 4 and Sr 2 Cu0 2 Cl 2 §, Ca 2 Cu0 2 Cl 2 §, 
NaV 2 5 |], Nd 2 Cu0 4 §, and 1-D insulators Sr 2 Cu0 3 
and Sr 2 Cu0 2 (|]. The measured signal is resonantly en- 
hanced by tuning the incident photon energy to lie near 
the Cu K or V L3 edge. The measurements have re- 
vealed remarkably similar characteristics as a function of 
photon energy loss: (1) the presence of a large, sharp 
and relatively dispersionless peak centered around a few 
eVs, and (2) the development of a low energy peak dis- 
persive towards higher frequencies for photon momentum 
transfers from the Brillouin zone (BZ) center along either 
the BZ edge or diagonal. Example of the data taken on 
Ca 2 Cu0 2 Cl 2 || is shown in Figure 1. The high energy 
peak has been associated with photon-induced charge 
transfer between orbitals of different atoms or differ- 
ent orbitals of the same atom[|| |6|, while the low fre- 
quency peak has been associated with a transition from 
the lower to upper Hubbard band across an effective Mott 
gap(|, |],J6| and a g-dependence of the Mott gap has been 
inferred However it does not seem obvious why an ex- 
citation across a Mott gap would show dispersion given 
that the physics of the Mott transition is highly local in 
character. 

Theoretical calculations on inelastic X-ray scattering 
have to our knowledge been limited to energy-band model 
calculations and exact diagonalization studies of small 
clusters [Q. While energy-band calculations might be ap- 
propriate for ground state properties of weakly corre- 
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FIG. 1: Experimental data for Ca2Cu02Ci2 obtained in 
Ref. U for momentum transfers along the BZ diagonal (a) 
and BZ edge (b), respectively. The values of the parameter 
X — [cos(q x a)+cos(q y a)]/2 indicate the values of the momen- 
tum transfer (q x ,q y ). Note that for panel (b) the momentum 
transfer runs from (0, 0) to (n, 0) only, X — 1 to 0. 



lated systems they do not adequately address the role 
of intra-atomic electron correlations which crucially af- 
fect properties of the excitation spectra of strongly cor- 
related systems. Exact diagonalization studies of small 
clusters [|| largely focuses on the energy separation be- 
tween the states excitable by the X-rays such as excitons, 
and suffers the limitation that the lineshape of the cal- 
culated spectra depend sensitively on cluster size due to 
finite size effects on electron dynamics. Thus it is clearly 
of interest for formulate a theory for inelastic X-ray scat- 
tering which does not suffer from cluster-size effects and 
is able to properly account for intra-atomic electron cor- 
relations. 

Two important features of the experimental data have 
yet to be clarified. First, the selection rules coming from 
the different orientations of the polarization directions of 
the incoming and outgoing photons as well as the direc- 



tion of their scattered momenta have not been used to 
determine the symmetry of the upper Hubbard band, for 
example. These selection rules have led to intense inves- 
tigation of the dynamics of electrons in strongly corre- 
lated systems like the high temperature superconductors 
to determine information about charge dynamics on re- 
gions of the BZ or the symmetry of the order parameter in 
the superconducting state ||. Secondly, the mechanisms 
and/or pathways of electronic relaxation revealed by in- 
elastic X-ray scattering have yet to be exploited. More 
precisely, besides the magnitude of the matrix elements 
coupling the conduction band to the excited states via 
the photon vector potential, the way in which the X-ray 
induced charge imbalance relaxes has remain largely un- 
explored. Inelastic X-ray scattering affords an open win- 
dow to investigate the symmetry and pathway of charge 
dynamics in strongly correlated systems. This is the topic 
of this paper. 

We choose to focus attention on non-resonant X-ray 
scattering, or scattering in which the frequency depen- 
dence of the incoming or outgoing photons can be individ- 
ually neglected and only the frequency shift 17 = LUi — uj s 
enters, where 0Ji iS denotes incident, scattered X-ray en- 
ergies, respectively. This certainly means that we have 
lost the ability to make quantitative predictions concern- 
ing the overall intensity of the scattering and we cannot 
comment on lineshape changes induced by varying the in- 
cident photon frequency. However our goal is to evaluate 
inelastic X-ray scattering in a model in which the cor- 
relations can be handled exactly - the Falicov-Kimball 
model in infinite dimensions - to determine which fea- 
tures emerge from the strong correlations. 

The Falicov-Kimball model, which has been used to 
describe a variety of phenomenon in binary alloys^), 
contains itinerant band electrons and localized electrons, 
in which the band electrons can hop with amplitude t* 
between nearest neighbors and interact via a screened 
Coulomb interaction U with the localized electrons: 



H = ~27<L ^ iCj+Ef ^ Wl ~^ ^ 4 c i+ U Yl c i c i w i> 

* * 1 

(1) 

where cj , Ci is the spinless conduction electron creation 
(annihilation) operator at site i and Wi = or 1 is a classi- 
cal variable of the localized electron number at site i. Ep 
and y, control the filling of the localized and conduction 
electrons, respectively. This is solved by dynamical mean 
field theory as described in detail in Ref. jl(| , where the 
reader is referred to for details. 

In this single band model with energy e(k), the in- 
elastic X-ray response is given formally by an effective 
density-density correlation function S(q, oS) = — i[l + 

™M]x"(q,^) with 

X (q, W ) = <[p(q),p(-q)]) H (2) 
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FIG. 2: Coupled Dyson equations for the inelastic X-ray 
scattering density-density correlation functions described by 
the scattering amplitude y a . Panel (a) depicts the Dyson 
equation for the interacting correlation function, while panel 
(b) is the supplemental equation needed to solve for the corre- 
lation function. The symbol F stands for the local dynamical 
irreducible charge vertex given in Eq. . In situations where 
there are no charge vertex corrections, the correlation func- 
tion is simply given by the first (bare-bubble) diagram on the 
right hand side of panel (a). 

formed with an "effective" density operator given by 

P(q) = J2 7a(k)4(k + q/2)c CT (k - q/2). (3) 

k,o- 

The strength of the scattering 7 a is determined by the 
curvature of the band as 

Here e 1,s denote the incident, scattered X-ray polar- 
ization vectors, respectively, and we have chosen units 
ks = c = h = 1 and have set the hypercubic lat- 
tice constant equal to 1. We can classify the scatter- 
ing amplitudes by point group symmetry operations. If 
we choose e l — (1,1,1,...) and e s = (1,-1,1,-1,...), 
then we have the Bi g sector, while e 1 = e s = (1, 1, 1, ...) 
projects out the A\ g sector since the i?2 3 component is 
identically zero in our model due to the inclusion of only 
nearest-neighbor hopping. We thus can cast the scatter- 
ing amplitudes in a simple form: 7^ (k) = — e(k) and 

iB ls (k) = t* E^cosk^-ip/VZ 

The Dyson equation for the density-density correlation 
function takes the form given in Fig. Note that there 
arc two coupled equations illustrated in Figs. || (a) and 
(b); these equations differ by the number of 7 a factors in 
them. The irreducible vertex function T is the dynamical 
charge vertex |nj] which takes the form 

1 (iuj m ,iuj n ;ivi^o) = 6 mn — — — . (5) 

J urn ~ Cjjti+; 

on the imaginary axis [iu) m = i7rT(2m + 1) is the 
Fermionic Matsubara frequency and iv\ = linTl is the 
Bosonic Matsubara frequency]. Here S m = ~E(iuj m ) is the 




i — i — i — ri — i — i — 
// X x=-1 




1234567012345678 

Frequency [t*] 

FIG. 3: Inelastic X-ray scattering response in the Bi g chan- 
nel along (a) the Brillouin zone diagonal and (b) along the 
zone edge for the half-filled Falicov-Kimball model on a hy- 
percubic lattice. The solid, dotted, short-dashed and long- 
dashed curves correspond to temperatures T = 1 .0, 0.5, 0.25, 
0.1, respectively. 
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FIG. 4: Inelastic X-ray scattering response in the A\ s chan- 
nel along (a) the Brillouin zone diagonal and (b) along the 
zone edge for the half-filled Falicov-Kimball model on a hyper- 
cubic lattice. The solid, dotted, short-dashed and long-dashed 
curves correspond to temperatures T — f.0, 0.5, 0.25, 0.1, re- 
spectively. 



local self energy on the imaginary axis and G m = G(iio m ) 
is the local Green's function on the imaginary axis, ff the 
scattering amplitude 7 does not have a projection onto 
the full symmetry of the lattice, then there are no vertex 
corrections from the local dynamical charge vertex jl2| . 

A straightforward calculation, shows that the B lg re- 
sponse has no vertex corrections on the zone diagonal 
q = (q,q,q,q, ...). Hence, the Bi g response is the bare 
bubble: 

XBi g (q,^) = ^ J du{f{u)xo{w,X,v)- f{uj + v) 

x Xofa X, v)[f{u) - f(uj + v)]xq(u}; X, v)} 

(6) 



with 

Xo(u;X,i>) = - j <!(/,{< ) 

x R 

and 

Xo(u;X, v) = - I <l<i>« ) 

x R 
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(8) 



Here we have used the following notation: f(w) = 
1/[1 + exp(w)] is the Fermi-Dirac distribution, p(e) — 
exp(— e 2 )/y / 7r is the noninteracting density of states; 

is the local self energy on the real axis; X — 
lim^co Yli cos qi/d = cos(q) for the zone-diagonal 
wavevector q = (q,q,q, ...,q); and F oa (z) = J dep(e)/(z — 



e) is the Hilbert transform of the noninteracting density 
of states. Techniques for finding the self energy |f(| have 
appeared elsewhere. 

The Ai g response everywhere and the B\ g response off 
of the zone diagonal, do have vertex corrections. The 
calculation of each response function is straightforward, 
but tedious. One needs to first solve the coupled equa- 
tions depicted in Fig. || on the imaginary axis and then 
perform the analytic continuation as in the Raman scat- 
tering case fL3|| . The end result is cumbersome and will 
be presented in a longer paper. 

The results for a correlated insulator U = 4t* at dif- 
ferent temperatures are shown in Figs. || and |I] for 
Big and Ai g inelastic X-ray scattering, respectively, as 
a function of transfered energy for different momentum 
transfers throughout the BZ measured by the factor X. 
Panel (a) for Figs. || and f| refer to scattering along the 
zone diagonal X = cos q for the zone-diagonal wavevec- 
tor q = (q,q,q, ...,q), and panel (b) refer to scattering 
along the zone edge [here we have q = (q, 0, q, 0, q, 0) 
for 1 > X = (1 + cos q) j '2 > and q = (-7T, q, ir, q, it, q) 
for > X = ( — 1 + cosg)/2 > — 1 The curves have been 
shifted vertically for clarity. The lowest set of curves 
X = 1 correspond to Raman scattering with optical pho- 
tons fL3|| , The main qualitative feature in both Figures 
are the presence of a small, dispersive low-energy peak 
for frequencies ~ t* and a large, dispersionless charge- 
transfer peak ~ U. While the charge-transfer peak re- 
mains relatively robust with increasing temperature, the 
low energy peak gains intensity from zero as temperature 
is increased. In particular all momenta show the devel- 
opment of low-energy spectral weight as T increases and 
there is a non-dispersive isosbestic point - a frequency at 
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FIG. 5: Plots of the low energy (panel a) and charge transfer 
(panel b) peak positions (solid symbols) and broadening (full 
width at half maximum, open symbols) for T = t* determined 
from Figs. 3 and 4 for Ai g , Bi g (circles, squares), respectively. 

which the spectra are temperature independent - around 
~ U/2. The high energy peak reflects the energy scale 
for excitations across the Mott gap and is relatively dis- 
persionless due to the local nature of the correlations. 
In contrast the low energy feature is a consequence ther- 
mally generated double occupancies which open up a low 
energy band up to energies ~ t* able to scatter X-rays. 
For decreasing temperature the low energy intensity dis- 
appears and only scattering across the Mott gap remains. 

The charge-transfer peak is sharp near the BZ cen- 
ter (X = 1), but broadens for momentum transfers ap- 
proaching the BZ corner X = — 1, more so for the B\ g 
channel than for A\ g . In fact A\ g and B\ g are identical 
at the (7r, 7r, 7r) point X = — 1 due to the local approx- 
imation. Any variation in the signal at the zone corner 
in different symmetry channels must be due to nonlocal 
many-body correlations. 

An important difference is that the Ai g results have 
no low-energy spectral weight for q = 0, corresponding 
to inelastic Raman scattering The vertex corrections 
remove all remnants of the low-energy response here, but 
it enters for any finite value of q. For an unpolarized 
measurement, the X-ray response is a superposition of 
the Bi g and A\ g spectra. 

We plot the behavior of the peak position and peak 
width (full-width at half maximum) for both the low- 
energy peak and the charge-transfer peak for both chan- 
nels in Fig. |5|. One can see more clearly that the low- 
energy peak has a width larger than its energy for both 
channels and for the A\ g channel follows the behavior of 
the corresponding B\ g feature away from the zone cen- 
ter. The charge transfer peak on the other hand is well 
defined for both channels. The only dispersive feature of 
the charge-transfer peak is the width of the B\ g peak. 

Referring back to the experimental data shown in Fig. 
El it is tempting to associate the relatively dispersion- 



less high energy peak with an excitation across a charge- 
transfer gap and the broad low energy peak with the 
dispersive feature generated from double occupancies. 
However the experimental data is not yet complete as 
the polarization and temperature dependence have not 
been measured. Our theory would predict, if this in- 
terpretation were correct, that the low energy feature 
would decrease in intensity as temperature is lowered. 
Moreover a polarization-dependent measurement could 
perhaps deconvolve the high-energy peak into two sep- 
arate peaks of A% g and B\ g symmetry, and would also 
be able to separate different behavior of the low energy 
peak near the zone center. Thus we believe it would be 
quite interesting to examine inelastic X-ray scattering at 
well-controlled temperatures and with polarizers for the 
incident and scattered light. We believe that a number 
of new and interesting features of charge excitations in 
correlated systems are likely to emerge if this can be ac- 
complished. 
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